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Abstract 

Cross sections for elastic production of J mesons in photoproduction and electropro- 
duction are measured in electron proton collisions at HERA using an integrated lumi- 
nosity of 55 pb~^. Results are presented for photon virtualities up to 80GeV^. The 
dependence on the photon-proton centre of mass energy VF^p is analysed in the range 
40 < Wyp < 305 GeV in photoproduction and 40 < W^p < 160 GeV in electroproduc- 
tion. The W^p dependences of the cross sections do not change significantly with and 
can be described by models based on perturbative QCD. Within such models, the data show 
a high sensitivity to the gluon density of the proton in the domain of low Bjorken x and low 
. Differential cross sections da/dt, where t is the squared four-momentum transfer at 
the proton vertex, are measured in the range |t| < 1.2 GeV^ as functions of W^p and Q^. 
Effective Pomeron trajectories are determined for photoproduction and electroproduction. 
The J/ t(j production and decay angular distributions are consistent with s-channel helicity 
conservation. The ratio of the cross sections for longitudinally and transversely polarised 
photons is measured as a function of and is found to be described by perturbative QCD 
based models. 
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1 Introduction 



Quantum Chromodynamics (QCD), the field theory of quark and gluon interactions, is expected 
to describe the strong force between hadrons. QCD is a successful theory in the limit of short 
distances, corresponding to small values of the strong coupling constant ag, where perturbative 
methods can be applied (perturbative QCD, pQCD). The bulk of the scattering cross section of 
hadrons however, is dominated by long-range forces ("soft interactions"), where a satisfactory 
understanding of QCD still remains a challenge. A large fraction of these soft interactions is 
mediated by vacuum quantum number exchange and is termed "diffractive". In hadronic inter- 
actions, diffraction is well described by Regge theory, where it is due to the t-channel exchange 
of a leading trajectory with vacuum quantum numbers, called the "Pomeron" trajectory. In the 
high energy limit, Pomeron exchange dominates over all other contributions to the scattering 
amplitude and leads to an almost energy-independent total cross section. Elastic photoproduc- 
tion of vector mesons, jp VM p, is a particular example for a diffractive process. Measure- 
ments of the cross sections for the elastic production of light vector mesons (p, a;, and 0) in low 
Q"^ electron-proton collisions at HERA as function of the photon-proton centre of mass energy 
W^p [1,2] have verified the expected universal Regge behaviour. 

The cross section for elastic photoproduction of J/tjj mesons, 7p — > J/ipp, on the contrary, 
rises steeply with Wryp [3-6], incompatible with a universal Pomeron. Due to the large mass 
of the J meson, which provides a "hard" scale (equivalent to a short range of the forces 
involved), the elastic photoproduction of J/^ mesons is expected to be described by pQCD. In 
electroproduction the photon virtuality can provide a second hard scale in addition to the 
J/ip mass, allowing the interplay between these two scales to be studied. 




Figure 1 : Elastic J/ ip production, a) in an approach based on Pomeron (P) exchange and b) in 
a pQCD approach via two gluon exchange. The kinematic variables are indicated in a). 

The elastic production of J/^ mesons is illustrated in figure 1. In QCD at lowest order the 
process is mediated by a colour-singlet state of two gluons (figure lb) and the cross section is 
related to the square of the gluon density in the proton. The gluon's momentum fraction x is 
kinematically related to W^p : the steep rise in the gluon density towards low values of x thus 
explains the steep rise of the cross section with increasing W-yp observed in the data. Beyond 
this approximation correlations between the gluons have to be taken into account and the cross 
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section for elastic J/ production involves the generalised gluon density (see for example [7] 
for a review). 

The dependence of the elastic J/ t/j cross sections on the squared four-momentum transfer t at 
the proton vertex shows a fast fall with increasing |t| . This dependence can be parameterised as 
an exponential function e^* at low values of although other shapes have also been proposed 
(for example [8]). In Regge theory with a t dependent Pomeron trajectory, the t dependence of 
the cross section varies with W^p, the slope parameter b increasing logarithmically with W^p 
("shrinkage" of the diffractive peak). In QCD-based models, on the other hand, the dependence 
of b on Wjp is expected to be weak [9, 10]. In addition to the elastic process, in which the 
proton remains intact, diffractive J/?/; production can lead to proton dissociation, 7p J /ipY , 
in which a low mass baryonic state Y is produced. This process is expected to be important at 
large values of\t\. 

In the past, diffractive J/ip cross sections have been measured using photon and electron beams 
in fixed target experiments up to centre of mass energies of about 20 GeV. At HERA the kine- 
matic range is extended up to photon-proton centre of mass energies of W^p ~ 290 GeV in 
photoproduction, and in electroproduction up to photon virtualities of < 100 GeV^ [1,3- 
6, 1 1-16]. In this paper new data are presented on the Q^, W-yp and t dependence of the cross 
section for elastic J/V' production. The data correspond to a factor of three more integrated 
luminosity than our previous publication for photoproduction [4] and a factor of two more for 
electroproduction [11]. The kinematic range is extended to values of W^p up to 305 GeV in 
photoproduction, while in electroproduction the range covered is 40 < W.yp < 160 GeV. Fur- 
thermore, the angular distributions for production and decay of the J / ijj mesons are determined 
in order to extract the cross sections of longitudinally and transversely polarised photons and to 
test the hypothesis of s-channel helicity conservation (SCHC), which predicts that the helicity 
of the J/ meson in the final state is the same as that of the initial (virtual) photon. 

2 Models for Elastic J/^ Production 

Within the Regge framework (see for example [17] for a review) the cross section for diffractive 
photoproduction of vector mesons at low values of |t | approximately follows a power law a^p oc 
W^p with 5 ^ 0.2 [18]. The power is related to the Pomeron trajectory 5 ^ A {aip{t) — 1) where 
c^ip (t) = ajp (0) + a'jpt. The existing measurements for J/ ip mesons, however, indicate a much 
steeper dependence on Wjp {5 ~ 0.8) than is predicted by the universal ("soft") Pomeron. There 
are also indications that the slope a'^p, responsible for the shrinkage of the diffractive peak of 
elastic J/'tjj photoproduction, is smaller [4] than the value of 0.25 GeV^^ expected from the 
soft Pomeron trajectory. To overcome this difficulty, Donnachie and Landshoff have suggested 
an additional "hard" Pomeron trajectory [19] for processes which involve a hard scale, such as 
the vector meson mass or a large momentum transfer Q^. With this conjecture the concept of 
a single universal Pomeron trajectory becomes obsolete for hard scattering processes. It has 
become customary, however, to introduce an effective Pomeron trajectory a{t) = + a't, 
where the intercept ao can be calculated within certain QCD models (see for example [20] for 
a review). 
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In photoproduction of J/ if} mesons the mass Mj/^ may serve as a hard scale and in electro- 
production both Mjj^ and Q"^. A third hard scale may be provided by a sufficiently large 
momentum transfer |t| at the proton vertex [13, 14]. In the presence of a hard scale QCD fac- 
torisation methods (e.g. coUinear factorisation, kx factorisation) may be applied. Factorisation 
allows the separation of the scattering amplitude into a perturbative hard scattering coefficient 
function and non-perturbative quantities, such as the input gluon density for the proton and the 
vector meson wave function. 

Early pQCD predictions, for example [21], assume that the two exchanged gluons have the same 
longitudinal momentum fraction x with respect to the proton, where x ~ (Q^ + M'^)/{Cp -\- 
W^p), and that each of the quarks making up the J/ijj meson carries half of the photon mo- 
mentum. Such models correspond to a leading approximation in log l/x, and at high energies 
the cross section depends on the square of the gluon density within the proton. More recently 
generalised or "skewed" parton distributions have been considered, where the two gluons have 
different fractional momenta [22-26]. 

The data presented here are compared with a pQCD model by Martin, Ryskin and Teubner 
(MRT [26]) which is based on kr factorisation and uses a parton-hadron duality ansatz avoid- 
ing the large uncertainties from the poorly known J/^p meson wave function. In this model, 
effects beyond the leading logarithmic approximation in log are included at the amplitude 
level, requiring an integration over the transverse momenta of the two gluons and hence the use 
of unintegrated gluon distributions. In the MRT calculations these distributions are derived [25] 
from the conventional integrated parton distributions \ as extracted from inclusive deep-inelastic 
scattering. The skewing effects are estimated independently by applying a factor to the ampli- 
tude [27]. Since these calculations only apply to the imaginary part of the scattering amplitude, 
dispersion relations are used to estimate the effects of the real part. In the parton-hadron duality 
approach the correct spin-parity state (J^ = 1 ) of the cc pair is projected out by using the 
appropriate rotation matrices in the integrals over the resonance mass region. Since the choice 
of the mass range in the integration is arbitrary to some extent, the normalisation of the cross 
sections is predicted with limited accuracy. The overall normalisation contains additional uncer- 
tainties due to missing higher order corrections. The approximations used are however believed 
to have little influence on the W^p and Q'^ dependences of the cross sections [28]. Predictions 
are provided both in the photoproduction and electroproduction regimes. 

The calculations by Frankfurt, McDermott and Strikman (FMS [9]) are based on the dipole 
approach. Here the exchanged photon turns into a qq pair long before the interaction with the 
proton. A leading logarithmic approximation for the interaction of this qq pair, described as a 
small transverse- size dipole, is used. For the interaction with the proton two-gluon exchange 
is assumed. In addition the effect of a running quark mass, a W^p dependent slope of the 
exponential t distribution, and generalised gluon distributions are considered in this calculation. 
Similarly to the MRT calculations the model does not provide an accurate normalisation of the 
cross section. Predictions are only available for photoproduction. 

'in this model the unintegrated gluon distribution is determined from the derivative of the standard gluon 
distribution with respect to log Q^, i.e. essentially from the second derivative of the proton structure function 
F2{x, Q^) with respect to log Q^, which in the kinematic region of the J/ip analysis, at low x and Q^, is not well 
measured at HERA yet. 



6 



3 Data Analysis 



The data were recorded with the HI detector in the years 1999 and 2000 when HERA was 
operated with electrons or predominantly with positrons^ of 27.5 GcV and protons of 920 GcV. 
The J /ip mesons are detected via their decays into /x+yU^ or e+e~ pairs. They are selected from 
data corresponding to an integrated luminosity of 55 pb^. 

3.1 The HI Detector 

The HI detector is described in detail in [29]. Charged particles are detected in the central and 

forward^ tracking detectors (CTD and FTD), which consist of drift and proportional chambers 
that provide a polar angle coverage between 7° and 165°. Tracks at large 6 are detected in the 
backward silicon tracker (BST [30], 165° < 9 < 175°). The central Uquid argon (LAr) [29] 
and backward lead scintillator (SpaCal) calorimeters [3 1 , 32] cover the polar angle regions 4° < 
e < 153° and 153° < 9 < 177.5°, respectively. For > 2 GeV^ the scattered positron is 
detected in the SpaCal, while the decay electrons from the J/ if) meson are identified in the 
LAr and SpaCal calorimeters. Muons are identified as minimum ionising particles in the LAr 
calorimeter or in the instrumented iron return yoke of the solenoidal magnet which surrounds 
the central detector (central muon detector, CMD, 4° < 6* < 171°). 

Dissociated proton states Y with masses My > 1.6 GeV may, after a secondary interaction, be 
measured in a set of detectors in the forward direction. These are the proton remnant tagger 
(PRT), an array of scintillators covering 0.06° < 9 < 0.17°, the drift chambers of the forward 
muon detector (FMD) [33] closest to the beam interaction region in the angular range 3° < ^ < 
17° and the forward region of the LAr calorimeter {9 < 10°). 

HI uses a multi-stage trigger system. At level 1 signals from the CTD, SpaCal, and CMD are 
used to obtain the present data sets. At level 2 information from these detectors and the LAr 
calorimeter is used in neural network algorithms [34]. 

The luminosity is determined from the rate of Bethe Heitler events. 

3.2 Kinematics 

The kinematics of the process ep epj /iIj are described by the following variables: the square 
of the ep centre of mass energy s = (p + A;)^; the negative four- momentum transfer squared 
at the lepton vertex = —q^ — —{k — k'Y; the four-momentum transfer squared at the 
proton vertex t — {p — p'Y and the inelasticity y — {p ■ q)/{p ■ k). The four-momenta k, k', 
p, p' and q refer to the incident and scattered positron, the incoming and outgoing proton (or 
dissociated system Y) and the exchanged photon, respectively. The centre of mass energy of 
the photon-proton system W.yp is given by W'^^ = {p + (if = ys — neglecting the proton 
mass. 

^Hereafter the term 'positron' is used for all lepton beam particles, whereas 'electron' is used for both electrons 
and positrons from J/ tp decays. 

^The positive 2;-axis is defined by the proton beam direction. The polar angle is measured with respect to the 
z axis and 6 < 90° is called the 'forward' direction. 
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In electroproduction the event kinematics are reconstructed using the double angle method [35], 

sin^^e(l — cos^^^) 
sin 6^ + sin Of. — sin(^^e + Sip) 
sin6'^(l + cos^e) 
^ sin 9^ + sin 9e — sin(6'e + 9^) 

Here is the energy of the incident positron and 9^, and 9^ are the polar angles of the J/ip 
meson and the scattered positron, respectively. The variable t is calculated as t ~ —{pt,ip+Pt,ey, 
where pt,^ is the transverse momentum of the J/ ijj meson candidate and pt^e that of the outgoing 
positron. 

In photoproduction, where the positron is not observed in the central detector, y is reconstructed 
via y — {E — pz)^/{2Ee) [36], where E and Pz denote the energy and the longitudinal com- 
ponent of the momentum of the J/ ^ meson. The variable t is approximated as i ~ ~Pt,ip (^^^ 
also the section on da/dt below). 

3.3 Monte Carlo Simulation 

Monte Carlo simulations are used to calculate acceptances and the efficiencies for triggering, 
track reconstruction, event selection and lepton identification. 

The elastic J/'tp signal events are generated using the program DIFFVM [37] which is based 
on the Vector Dominance Model and permits separate variation of the dependence on W^p, t 
and . The parameters are iteratively adjusted to those of the present measurements. DIF- 
FVM is also used to generate J/ tp production with proton dissociation. A mass dependence of 
da/dM^ oc /(M|-) My ^ is implemented, where /(M|-) = 1 for > 3.6 GeV^. At lower 
My the function / {My ) takes into account the production of excited nucleon states. The decay 
angular distributions of the J/tjj meson are simulated assuming s-channel helicity conservation. 
For electroproduction, radiative corrections are included using the generator HERACLES [38], 
where contributions up to order aq^D taken into account. 

The non-resonant background is estimated using the generators LPAIR [39], which simulates 
the process 77 — > i'^l" and COMPTON [40] for the QED Compton process ep — > e^p. Cross 
checks with the generator GRAPE [41] did not show significant deviations from the results of 
LPAIR in the region of the present analysis. 

For all processes, detector effects are simulated in detail with the GEANT program [42]. The 
detector response including trigger efficiencies is tuned using independent data. Remaining 
differences are included in the systematic errors. The simulated events are passed through the 
same reconstruction software as the data. 

3.4 Event Selection 

Elastic J/ t/j events are selected by requiring two muons or two electrons and, in the case of 
electroproduction, a scattered positron candidate. For photoproduction the absence of any such 
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candidate is required. As described in table 1 four data sets are defined covering different 
regions of and W^p and corresponding to different signatures of the J/V' decay leptons. 

For data set I (electroproduction) the scattered positron must be detected with an energy of 
at least 12 GeV in the SpaCal and the reconstructed value of must be within 2 < < 
80 GcV^. To suppress photoproduction background and to reduce the fraction of events with 
initial state QED radiation, events are rejected if Y^{E — p^) < 45 GeV, where the sum runs 
over all final state particles including the scattered positron. Neglecting radiative effects this 
variable is expected to be twice the incident positron energy due to longitudinal momentum 
conservation. 

For the selection of photoproduction events (data sets II-IV) the absence of any candidate for 
the scattered positron is required, restricting the accepted range of negative four-momentum 
transfer squared to below about 1 GeV^, with (Q^) = 0.05 GeV^. 

Data sets I and n (40 < W^p < 160 GeV) contain J/^ — > events. Exactly two oppositely 
charged particles must be present in the CTD, with transverse momenta (with respect to the 
beamline) pt > 0.8 GeV. A reconstructed vertex within ±40 cm of the z coordinate of the 
nominal beam interaction point is required. At least one particle must be identified as a muon 
in the central calorimeter or in the CMD. For data set n background from cosmic ray muons is 
rejected using an acoUinearity cut as well as timing information from the CTD. Further details 
of this analysis may be found in [43]. 

Data sets III and IV are selected to cover photoproduction at high values of W^p, which are 
related to large polar angles of the J ftp decay leptons. The J/ijj decay into e+e" pairs is 
used. Data set 111 (135 < W^p < 235 GeV) requires one decay electron to be measured in 
the CTD coming from within ±40 cm of the nominal beam interaction point and one in the 



Data set 


I 


II 


III 


IV 


Kinematic region 


Electroproduction 


Photoproduction 


g2 range [GeV^] 


2-80 


< 1 


{Q')[GeV'] 


8.9 


0.05 


W^p [GeV] 


40 - 160 


135 - 235 


205 - 305 


\t\ [GeV'] 


< 1.2 


Decay channel 




J/ip ^ e+e 


Lepton signature 


Track-Track 


Track-Cluster 


Cluster-Cluster 


Lepton polar 
angle region [°] 


20 - 160 


^1 : 80 - 155 
e-2 : 160 - 177 


01 : 160 - 174 
62 : 160 - 175.5 


Lepton energy 

[GeV] 


Pt > 0.8 


Pt,i>0.7,pi >0.8 
E2 > 4.2 


Ei^2 > 4.2 
max(£'i, E2) > 6 


Elastic selection 


No signal in forward detectors 


jLdt [pb-^] 


54.79 


30.26 


26.90 



Table 1: Summary of the most important event selection criteria for the four different data sets 
together with the corresponding integrated luminosities. 
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backward calorimeter SpaCal. The selected polar angle regions are given in table 1. The elec- 
tron measured in the CTD must have a momentum pi > 0.8 GeV and a transverse momentum 
Pi 1 > 0.7 GeV and must be identified by a matching electromagnetic energy deposition in the 
central calorimeter. The other electron is selected by requiring a cluster in the SpaCal with an 
energy E2 > 4.2 GeV. Data set IV (185 < Wyp < 305 GeV) requires both electrons to be 
detected as energy clusters in the SpaCal with energies £'12 > 4.2 GeV and the more energetic 
cluster to be above 6 GeV. At least one electron must be in the acceptance region of the BST 
and every electron in the BST acceptance region must be validated by a BST track from the 
nominal interaction point. This requirement rejects most of the non-resonant background from 
Compton scattering. In both data sets in and IV the energy in the SpaCal outside the selected 
electron cluster(s) must be negligible. Further details of this analysis may be found in [44]. 

In order to suppress background from proton dissociative or inelastic J/ production, no ad- 
ditional tracks are allowed in the CTD or FTD and the selected events are required to have no 
significant signals in the forward detectors (PRT, FMD and LAr). The fraction of proton disso- 
ciation is further suppressed by limiting t to the range \t\ < 1.2 GeV^, where elastic processes 
are dominant. These requirements reject most of the proton dissociative background. The re- 
maining fraction is 14% on average, ranging from 8% at \t\ to 35% at \t\ ^ 1.2 GeV^. 
It is corrected for using the MC simulation, which is tuned to give a good description of the 
forward detectors. A further correction is applied to account for ^(25') decays into J/ if) and 
neutral mesons. This correction is estimated to be 4% for data sets 1 and 11 and approximately 
2% in sets III and IV, where the neutral decays are partly rejected by the cut on the energy in 
the SpaCal. 

Triggers based on muon and track signatures from the decay leptons are used for data sets 1 and 
11. For data set 1 a trigger signal is also derived from the scattered positron. The triggers for 
data sets HI and IV are based on signals due to the J/ '4> decay electrons from the SpaCal and 
the CTD (set III) . In addition the triggers for data sets II-IV use second level triggers based on 
neural network algorithms. 

Figure 2 shows the two-lepton invariant mass distributions for the four data sets. The shapes 
of the J/ signal peaks reflect the usage of different detectors with different resolutions and 
a different response to electrons, muons and photons. The signal of data set III shows a tail 
towards low masses due to radiative energy losses of the electron reconstructed in the tracking 
detector. In all data sets the non resonant background below the J/ ip signal peaks is dominated 
by 77 t^l", where one photon originates from each of the positron and the proton. At 
high Wryp a potential source of background is Compton scattering ep e'yp where the final 
state electron and photon can form an invariant mass of the same order as the J/ip mass. It 
is efficiently suppressed by the BST track requirements explained above. In addition the BST 
tracks lead to an improved mass resolution in data set IV. 

The number of J/ ip events is determined in each analysis bin by a fit of the sum of a signal 
and a background function to the dilepton mass distribution. For data sets I and II {J/ip 
jJi^jj.') the signal shape is a Gaussian function, and the background is fitted using a power 
law distribution. For the signal in data set III (J/^ e~^e~) the radiative tail is taken into 
account by fitting a modified Gaussian distribution, f(Mee) oc ^ exp (— (Mee — ^)'^/{2a')) 
where a' = {a + r {\Mee — /i| — Mee + lAY- Here fi and a denote the peak position and 
the standard deviation and r parameterises the contribution of the radiative tail. In data set IV 
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( J/V' — > e+e") a single Gaussian function is adequate to describe the signal. In data sets in and 

IV the shape of the background is found to depend strongly on W^p. The shapes are reasonably 
well described by the predicted shapes of the Monte Carlo simulations LPAIR and COMPTON, 
which are therefore used in the fit. 

Data and MC simulation are compared in figure 3. Each row corresponds to one of the four 
data sets. The selected events from a mass window around the nominal J/ijj mass (±0.2 GeV 
in data sets I and II, ±0.3 GeV in data set IV and 2.6 < M^e < 3.4 GeV for data set IE) 
are shown before applying the cuts on the forward detectors. For data sets I and n the non- 
resonant background, which is small in this kinematic region, has been subtracted and the data 
are described by a combination of DIFFVM for elastic and proton dissociative J/ ijj production. 
For data sets III and IV the selected events are shown but without subtracting the non resonant 
background. The data are seen to be reasonably well described by a sum of simulations for 
elastic and proton dissociative J/ V' production (DIFFVM), 77 — > e+e~ (LPAIR) and ep — > e7p 
(COMPTON). 

3.5 Systematic Uncertainties 

The systematic uncertainties on the cross sections are dominated by detector effects which 
are not perfectly modelled in the Monte Carlo simulation. Most uncertainties are obtained 
by comparisons of data with simulation after tuning the detector simulation with independent 
data sets. The uncertainties on the measured cross sections are then estimated by variations 
of the simulation. In the following the main sources of the uncertainties are summarised and 
typical values are given for the uncertainty on the total cross section. 

• The uncertainty due to the track reconstruction efficiency in the CTD is 1% per track. 
The track information from the BST has two sources of uncertainty: coherent signal 
losses (3.0%) and track reconstruction efficiency (1.5%). 

• The uncertainty on the lepton identification efficiency leads to a cross section uncertainty 
of 1.5% for muons and 2% for electrons measured in the CTD. The uncertainty on the 
energy measurement of the decay electrons in the backward calorimeter is estimated to 
vary linearly from 2.7% at 3 GeV to 0.5% at 27.5 GeV from an analysis of Compton 
scattering [45]. The resulting uncertainties on the cross sections vary from 1% to 7%, 
depending on W^p. A small additional uncertainty for data set IV arises due to an uncer- 
tainty of 0.3 mrad in the reconstruction of the polar angle of the decay electrons in the 
BST, leading to a W^p dependent cross section uncertainty of 1 — 3%. 

• The uncertainties of the trigger efficiencies are determined to be 1.6%, 5%, 6.5% and 5% 
for data sets 1 to IV, respectively. 

• The separation of elastic events from proton dissociation leads to a systematic uncertainty 
of 4 — 6% due to the modelling of the response of the forward detectors, with a small 
dependence on W^p and \t\. The error due to the simulation of the dependence of the 
cross section on My was found to be negligible by comparison. 
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• The uncertainty in the modeUing of the z position of the interaction region affects the 
Wyp dependence of the cross section and is found to be 1% on average for data sets I and 
II, 0.5 - 2.6% for III and 2.0% for IV. 

• Varying the methods of determination of the number of signal events (e.g. by using a 
counting method instead of fits, or by changing the shapes of the background functions), 
results in a 1% uncertainty for data sets I and II (fi^ii^). For data sets III and IV (elec- 
trons) an uncertainty between 3% and 6% is estimated, which is due to the uncertainties 
in the signal and background shapes. 

• For the electroproduction sample, an additional uncertainty of 4% is estimated which 
covers uncertainties in the reconstruction of the energy and angle of the scattered positron. 

• Other sources of systematic uncertainties are the luminosity measurement (1.5%), the 
J/ijj branching ratio (1.7%) and the '4>{2S) background (0.5% for sets I and II, 1.5% on 
average for IE and IV). 



The systematic uncertainties are calculated in each analysis bin and the total uncertainty is 
obtained by adding all individual contributions in quadrature. The average values for the total 
systematic uncertainties on the cross sections are 8%, 9%, 10% and 11% for the data sets I to 
IV, respectively. The correlated part of the error, which affects all bins equally, is estimated to 
be approximately 5% and is not included in subsequent fits unless mentioned otherwise. 



4 Results 

Cross sections are calculated for the individual data sets I-IV using the number of selected 
events after correcting for non resonant, proton dissociative and ■'ip{2S) backgrounds as de- 
scribed in the previous section. The efficiencies A for the event selection are in general deter- 
mined from the MC simulation. In the equivalent photon approximation the 'yp cross section is 
given by: 

a(7P - J/#) = (1) 

Here $^ [46] denotes the photon flux in the and W^p range considered, C the integrated 
luminosity and BR the branching ratio for the decay of the J/t/j mesons'*. 

Note that this cross section corresponds to a^p = a'^p + ea^p, where cr^p and a^p are the cross 
sections for transversely and longitudinally polarised photons, respectively, and e is the po- 
larisation parameter of the virtual photon^. The parameter e depends only on the kinematics, 

e = {\ — y)/{\ — y + \y^)- In the kinematic range of the present analysis e is generally above 
0.95 with (e) = 0.993. Cross sections are given at 'bin centres', (W^7p), (Q^) ^nd (t), which are 
determined taking into account the measured W^p , Q"^ and t dependences. 

^Branching fractions (5.88 ± 0.10)% and (5.93 ± 0.10)% [47J are used for J/V' and e+e^, respec- 

tively. 

^The present results will be compared with results from the ZEUS collaboration [16], where g^.^ + al^^ is 
extracted. In the present kinematic region the difference is however small compared with the measurement errors. 
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4.1 Dependence 

The cross sections for elastic J/ ip production as a function of at W^p — 90 GeV are listed 
in table 2 and shown in figure 4a. The photoproduction point is obtained from the fit described 
in the next section. 

A phenomenological fit of the form u^p oc (M| + to the HI data yields a value of 

n — 2.486 ± 0.080(stat.) ± 0.068(syst.). This result confirms, with smaller errors, the 
dependence observed previously by HI [11]. The quality of the fit is good (x;^/ndf= 0.5). 
Recent results from the ZEUS collaboration [6, 16] are also shown in figure 4a, which agree 
well with the present data in the entire range of Q^. 

In figure 4b the pQCD calculations 'MRT' of Martin et al. [26] are compared with the fit result 
quoted above. Results with four different gluon distributions (CTEQ6M [48], MRST02 [49], 
HIQCDFIT [50] and ZEUS-S [51]) derived from global fits to current inclusive F2 measure- 
ments and other data are shown. A normalisation factor is determined individually for each 
prediction by comparing with the data across the complete range. The different factors, 
which are mainly given by the photoproduction measurement, are between 1.5 and 2.8. The 
theoretical predictions of the shape of the Q'^ dependence are consistent with the fit to the data 
within the experimental uncertainties, which are shown as a grey band in figure 4b. 

4.2 W^p Dependence 

The 7p cross section for elastic J / ijj production is presented as a function of W^p in figures 5a 
and 6a and in tables 3 and 4 for photoproduction and electroproduction, respectively. 

In figure 5a the photoproduction data are shown with the result of a fit of the form cr^p oc W^^. 
Separate relative normalisation factors for the three data sets are additional fit parameters which 
take into account the correlated systematic uncertainties. The fit yields a value of 5 = 0.75 ± 
0.03 ± 0.03. The first error is obtained using only the statistical uncertainties in the fit while the 
second one reflects the systematic uncertainties. The fit result is in agreement with our previous 
result [4]. Similar data from the ZEUS collaboration [6] (also shown in figure 5a) agree well 
with the present data. 

A comparison with theoretical predictions is shown in figure 5b, where the ratio of theory to 
the fit result is shown. The uncertainty of the fit result is indicated by the grey band. The MRT 
predictions are normalised using the factors obtained from the distributions. The same four 
gluon distributions are used to calculate the respective unintegrated skewed gluon distributions 
which are required by MRT. The W^p dependence is observed to be quite sensitive to the shape 
of the gluon distribution^. While the results based on the gluon distributions CTEQ6M and 
ZEUS-S describe the shape of the data well, the gluon distribution from the HI fit to inclusive 
data leads to a steeper Wyp dependence and the one from MRST02 to a flatter W^p dependence 
than is observed. The dipole model result EMS [9] based on the CTEQ4L [52] gluon density 
is somewhat too steep. Note, however, that these observations are based on the central val- 
ues of the respective gluon distributions and do not take into account their uncertainties. The 
kinematic range used in the MRT calculations extends to lower values of Bjorken x and Q'^ 

^For a detailed discussion of the sensitivity and the uncertainties of the model assumptions see [28]. 
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than was available in the inclusive data used for the determination of the gluon densities and an 
extrapolation to very low values of is performed. 

In figure 6a the electroproduction cross section is shown in three bins of (2 < < 5 GeV^, 
5 < < 10 GeV^ and 10 < < gO GcV^). Data from the ZEUS experiment [16], which 
are shifted to the present bin centres using the dependence measured by ZEUS, are in agree- 
ment. In figure 6a the results from MRT based on the gluon density CTEQ6M using the same 
normalisation factor as above are also shown and give a reasonable description of the data. 

The W^p dependence is found to be similar to that obtained in photoproduction. When param- 
eterised in the form W^^, the fits to the HI electroproduction data yield 5 values which are 
compatible with photoproduction within the rather large experimental errors (see table 5). The 

fitted values for 5 describing the IV^p dependence of elastic J /if) production from this analysis 
and from [6, 16] are displayed in figure 6b as a function of Q^. Within the present experimental 
accuracy no dependence on Q'^ is observed. 

4.3 Differential Cross Sections dcr /dt 

The t dependence of the elastic cross section for J / ijj meson production is studied in the range 
40 < W^p < 160 GcV for different bins. The differential cross sections da /At as derived 
from data sets 1 and 11 are listed in table 6 and shown in figure 7a with fits of the form da /dt oc 
e^^. The resulting 6 values (table 5) for electroproduction are systematically lower than the value 
for photoproduction but are compatible within the errors. 

In the context of developing the calculations using generalised parton densities, Frankfurt and 
Strikman [8] have proposed an alternative t dependence. It is based on a dipole function with 

a t dependent two-gluon form factor, leading to dcr/dt oc (1 — t/m|g)~^. In a fit to the pho- 
toproduction data the two-gluon invariant mass is left as a free parameter. A value of 
m2g = (0.679 ±0.006 ±0.011) GcV is obtained with xVndf = 5.5 compared to xVndf= 0.25 
for the exponential function. The dipole form is thus strongly disfavoured by the data. 

For photoproduction, the measurement of the t dependence has been extended to significantly 
higher W^p than in our previous publication [4] using data sets 111 (135 < W^p < 235 GeV) and 
IV (205 < Wjp < 305 GeV). Due to the reconstruction of the J/ t/j electrons via calorimeter 
signals the resolution in pf^, which is used to approximate t, is worse than in the track based 
measurements. The differential cross sections da/dt are obtained using an unfolding proce- 
dure [53]. The results (last two lines in table 7) are shown in figures 7b and c with exponential 
fits, which describe the data well. The resulting b values are listed in table 10 and are discussed 
further in the following section. 

4.4 Effective Pomeron Trajectories 

In models based on Regge phenomenology and Pomeron exchange, the energy dependence of 
the elastic cross section follows a power law: 
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where a{t) — ao + a't describes the exchanged trajectory and bo and Wq are constants. Equa- 
tion 2 relates the dependence of the differential cross section on t to that on Wjp by 



— (t) OC e^bo+4a'ln{W^j,/Wo))t_ 

Here only t dependent terms are kept. In hard interactions, where Regge phenomenology with a 
single universal Pomeron may no longer be applicable, an 'effective Pomeron trajectory' [20] is 
nevertheless often extracted in order to describe the dependence of the differential cross sections 
on Wyp and t. For the determination of this effective trajectory, a double differential analysis 
is performed in which the differential cross section da/dt is measured in bins of W^p and t. 
The measurements are displayed in figures 8a and b for photoproduction and electroproduction, 
respectively (tables 7 and 8). First, one-dimensional fits of the form W^p"*^^*^^"^^ to the cross 
sections in each | (t) \ bin are performed. The results, which are listed in table 9 and displayed as 
solid lines in figures 8a and b, describe the data well. In figures 9b and c the one-dimensional 
fit results for a{t) are compared with recent results [6, 16] from the ZEUS collaboration, which 
are in good agreement. 

A two-dimensional fit of the function given in equation 2 to the data yields values for bo, ao 
and a'. The parameter Wq is arbitrarily chosen to be 90 GeV; the fit result does not depend 
on this choice. As described before, different normalisations are allowed for the different data 
sets in the fit. Figure 9a shows the result of the two-dimensional fit for a{t) as solid and dashed 
lines for photoproduction and electroproduction, respectively. Error bands corresponding to one 
standard deviation are shown, taking the correlation between ao and a' into account. The results 
for a{{t)) from the one-dimensional fits are shown as points with error bars for comparison. 
Good internal consistency is observed. 

The results of the two-dimensional fits are listed in the following table. 



[GeV^] 


bo [ GeV-^] 


ao 


a' [ GeV-2] 


<1 
2-80 


4.630 ± 0.0601^63 
3.86 ±0.13 ±0.31 


1.224 ±0.010 ±0.012 
1.183 ±0.054 ±0.030 


0.164 ±0.028 ±0.030 
0.019 ±0.139 ±0.076 



Here the first errors are statistical and the second reflect the systematic uncertainties. 

The Wyp dependence of the cross section is predominantly determined by ao and the fit values 
lead to a W^p dependence very similar to the parameterisation with 6 discussed above. The 
parameter a' relates the t and W^p dependences and if non-zero leads to the 'shrinkage' of the 
diffractive peak. For photoproduction a' is larger than zero by four standard deviations and is 
two standard deviations below the value of 0.25 GcV~^ obtained for the soft Pomeron in [54]. 
For electroproduction a' is compatible with 0, which matches the expectation in [20], but due 
to the errors a' is also compatible with the value measured for photoproduction. 

Alternatively the value of a' can be measured using the dependence of the t slope parameter on 
Wjp , using b{W^p) = bo + Aa' ln.{Wjp/Wo). Exponential fits of the form e^^ to the measured 
differential cross sections da/dt in bins of W^p are performed and the resulting values for b are 
displayed in figure 10a and b and listed in table 10 for photoproduction and electroproduction. 
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For photoproduction the b values are seen to increase with W^p. These b values are independent 

of normalisation uncertainties between data sets. The curves in figure 10a and b show the 
corresponding result b{Wjp) from the two-dimensional fit described above. 

In figure 10a photoproduction results for the slope parameter from the ZEUS experiment [6] in 
a similar kinematic region are also shown. They show a similar dependence on W^p but are on 
average 0.5 GeV ~^ lower. This difference in the absolute size of b may be due to differences in 
the handling of the background from proton dissociative events, which has a much shallower b 
slope than for the elastic case (1.6 GeV~^ [13]). 

4.5 Helicity Studies 

The assumption that the J/^ meson observed in the final state keeps the helicity of the photon 
is referred to as s-channel helicity conservation (SCHC). This assumption can be tested by 
measurements of the angles in the production and decay of the J/ijj meson. If SCHC holds, the 
angular analysis leads to a separation of the cross sections due to longitudinally and transversely 
polarised photons which are both predicted in the MKT calculations. 



Three angles are defined, which are illustrated in the figure above. 6* is the polar angle of the 
decay muon with the charge of the beam lepton in the J/ ijj rest frame. 9* — 0° corresponds to 
the flight direction of the J/ t/j in the 7p centre of mass frame. 0* is the angle between the J/ t/j 
production plane, defined by the exchanged photon and the J/ ijj meson, and the decay plane 
in the centre of mass frame. $ is the angle between the scattering plane of the beam lepton 
and the J/ ip production plane. The angle $ can only be measured when the scattered electron 
is observed, i.e. for electroproduction. In the case of SCHC and natural parity exchange, the 
angular distributions of the J/ '4> production and decay are functions of cos 9* and ^ = 0* — $ 
only [55]. 

The angular distributions are expected to be similar for elastic and proton dissociative pro- 
cesses. This expectation is verified within the present statistical accuracy. Therefore, in order 
to increase the statistics, the cross sections differential in the angles are derived using a data 
set, which includes proton dissociative events in addition to data sets I and II. The selection of 
proton dissociative events is similar to that for data sets I and II (section 3.4 and table 1), but 
now a signal in one of the forward detectors is required and one additional track is allowed with 



yp centre of mass frame 



J/\|/ rest frame 




; J/\|/- production ^.^^^'^ \ 
plane 

J/\|f - decay 
plane 
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6 < 20°. Non-diffractive events are rejected by requiring z > 0.95, where z — Ej/^/Ej in the 
proton rest frame. Furthermore, an increased t range, |t| < 5 GeV^, is allowed, since the proton 
dissociative J/^p cross section shows a flatter t dependence, with an exponential slope of about 
1.6 GcV"^ [13]. 

In the present analysis, differential cross sections for the four angles are used to measure com- 
binations of seven of the 15 spin-density matrix elements, which describe the spin structure of 
the interaction completely^. The measured angular distributions and their dependence on the 
spin-density matrix elements r^xMXi^) i^^V- 



l+r°o'+(l-3rSo')cos^r (4) 
1 + r?ii cos(20*) (5) 
1 - £ri_i cos(2^) (6) 
1 - £ ('^00 + 2rn) cos 2$ + y/2e{l + e) (r^o + 2rn) cos (7) 
Here s is the polarisation parameter of the virtual photon. 

Figures 11a and b show the differential cross sections dcr/d cos 9* and da/dcj)* in four bins 
of Q^. Figures 11c and d show the differential cross sections da/d"^ and da/d<l> in three bins 
of Q^. The results of fits of equations 4, 5, 6 and 7 are shown as full lines. In the fits, the 
spin-density matrix elements or the combinations r^Q + 2rli and Tqc + 2rfi for equation 7 are 
free parameters. In figure 1 lb and d results from a fit assuming SCHC are also shown. 

The spin density matrix elements, which are determined by the fits, are shown in figures 12a- 
e as functions of Q'^ (tables 11 and 12). The analysis is also performed in bins of \t\ and 
the resulting spin density matrix elements are displayed in figures 12f-j (tables 13 and 14) as 
functions of 

Results from the ZEUS experiment [5,6, 15, 16] are also shown in figures 12a, b, c and f, which 
are in good agreement with the present results. In figures 12b, d, e, g, i and j the expectation 
from SCHC, namely 0, matches the data well. SCHC yields a relation between two spin density 
matrix elements: rl_^ = (1 — ^oo)/2- This is observed to be fulfilled within errors. 

In the case of SCHC, the matrix element Tqq provides a direct measurement of R, the ratio of 
the cross sections for longitudinal and transverse polarised photons, and a^, respectively: 



da 



d cos 9* 
da 

d0* 
da 

d^ 
da_ 

d$ 



oc 



oc 



oc 



oc 



^04 
00 



el 



00 



The values of R are presented in figure 13a and in table 11. For comparison the prediction 
from MRT [26] is shown, which depends only weakly on the gluon density. In figure 13a the 

^Spin density matrix elements ?'^(^);^(^) or ''5i'(7)a(i/j) linear combinations of the transition amplitudes 
'^xMMih) from a photon of helicity A(7) to a J/ip of helicity X{ip). 
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gluon density from CTEQ6M is chosen with the normaUsation as before, which gives the best 

description of the and Wyp dependences of the cross sections. The prediction is somewhat 
above the data but still describes the dependence reasonably well. Similar results from [6, 16] 
agree also with the present data. 

The values of R can be used to derive the cross sections and using the relationship a^p — 
cr^p + £cr^p. The results are shown in figure 13b as a function of Q"^. dominates at low Q^, 
while at ~ -M| both and are of similar magnitude. The MKT predictions are compared 
with the data using different gluon density parameterisations. The differences between the 
predictions are not very large. All gluon density parameterisations give a reasonable description 
of the data, although is somewhat above the data for > 3 GeV^. 

In brief, the helicity studies show consistency with SCHC within experimental errors. The ratio 
of cross sections for longitudinally and transversely polarised photons is extracted and its Q'^ 
dependence is found to be reasonably described by the MRT calculations. 

5 Summary 

New measurements are presented of elastic J/ ^ photoproduction and electroproduction in the 
ranges 40 < W^p < 305 GeV and 40 < W^p < 160 GeV, respectively \ 

The cross section a{'yp — > J/t/jp) is measured as a function of Q"^ in the range < < 
SOGcV^andafitoftheformcT^p oc (M| + Q2)-« yields a value of n = 2.486 ± 0.080(stat.) ± 
0.068(syst.). The shape of the distribution is well described by a perturbative QCD cal- 
culation by Martin, Ryskin and Teubner (MRT), almost independently of the gluon density 
distribution used. 

The photoproduction cross section is measured as a function of the photon-proton centre of 
mass energy W^p in the range 40 < Wjp < 305 GeV, and can be parameterised as a^p oc W^p 
with S — 0.754 ± 0.033(stat.) ± 0.032(syst.). The results for S in electroproduction, measured 
in the range 40 < W^p < 160 GeV, are consistent with those in photoproduction and no 
dependence is observed within experimental errors. Predictions of the W^p dependence of the 
cross section in pQCD-based models depend strongly on the gluon distribution, as can be seen 
explicitly in the MRT model. A good description of the shape of the data can currently be 
achieved only with some gluon parameterisations. This demonstrates the potential to constrain 
the gluon distribution with the elastic J/ip data in a kinematic region (low x, low ) where 
fits from inclusive data yield gluon distributions with large uncertainties. 

The differential cross section da/dt for elastic J/ip photoproduction for l^l < 1.2 GeV^ is mea- 
sured in the extended range of 40 < W^p < 305 GeV. A single exponential function yields 
a good description of da/dt in this range, while a functional form based on a dipole function 
is strongly disfavoured. The slope parameter b of the exponential shows a dependence on Wyp 
which is weaker than expected from soft Pomeron phenomenology, but is clearly positive, lead- 
ing to shrinkage of the diffractive peak. The slope parameter b in electroproduction agrees with 
the photoproduction values within errors, but has a tendency to decrease with increasing Q'^ . 

**The results of the present analysis agree within errors with our previous results. We consider the new data to 
supersede them due to improved statistics and better understanding of the detector efficiencies. 
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Effective Pomeron trajectories ckq + a't for elastic J /if) photoproduction and electroproduction 

are determined from a simultaneous analysis of dcr/dt as a function of W^/p and \t\. The elec- 
troproduction and photoproduction results are consistent with each other within errors. The 
trajectory for photoproduction has a t slope which is two standard deviations below the soft 
Pomeron value but four standard deviations above zero. 

Finally, the helicity structure of diffractive J/ip production is analysed as a function of Q"^ and 
\t\. No evidence is found for a violation of s-channel helicity conservation (SCHC). Assuming 
SCHC, the ratio of the longitudinal to the transverse polarised photon cross sections is deter- 
mined as a function of Q"^ and is found to be consistent with QCD calculations. 
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Figure 2: The dilepton invariant mass distributions (data and fits) in the four kinematic regions 
defined in table 1. 
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Figure 3: Observed event distributions for the four data sets (points) defined in table 1, omit- 
ting the forward detector cuts against proton dissociative events. The first two rows correspond 
to the selected J/i}) ^ f^^f^' candidates where the small non-resonant background has been 
subtracted. The data are shown with the elastic (signal) simulations (DIFFVM el., white area) 
and proton dissociation MC (DIFFVM pdiss., shaded area). Rows three and four correspond to 
J/ip ^ e~^e~ candidates, where the non-resonant background is not subtracted. Here, in addi- 
tion to the elastic and proton dissociative J/^ simulations the contributions from 77 e~^e~ 
(LPAIR) and Compton scattering (COMPTON) are shown. The normalisations are obtained 
from a fit of the overall mass peak of each data set. The variables W^p , Q^, t and p^^ are 
defined in the text. 9^ refers to the decay muons of data set II. In row three 9i and 62 refer to 
the decay electrons which are selected in different polar angular regions. In row four Ei and E2 
refer to the energies of the decay electrons. 
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Figure 4: a) Total cross section for elastic J/ijj production as a function of Q'^ in the range 
\t\ < 1.2 GeV^ at W^^p = 90 GeV. The inner error bars show the statistical errors, while 
the outer error bars show the statistical and systematic uncertainties added in quadrature. The 
solid line is a fit to the HI data of the form a^p oc (M^ + Q^y. Data from the ZEUS 
experiment [6, 16] are also shown, b) The ratio of the MRT calculations [26] to the fit from 
a). The MRT QCD predictions are based on different gluon distributions [48-51]. The curves 
are individually normalised to the measurements across the complete range yielding factors 
between 1.5 and 2.8. The shaded band represents the uncertainty of the fit result. 
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Figure 5: a) Total cross sections for elastic J/V" production as a function of W^y in the range 
|t| < 1.2 GeV^ in photoproduction. The inner error bars show the statistical errors, while the 
outer error bars show the statistical and systematic uncertainties added in quadrature. The solid 
line shows a fit to the HI data of the form a oc W^p. Results from the ZEUS experiment [6] 
in a similar kinematic range are also shown, b) The ratio of theoretical predictions to the fit 
to the HI data in a). The shaded band represents the uncertainties of the fit result. Predictions 
from MRT QCD calculations [26] and a dipole model (FMS, [9]) based on different gluon 
distributions [48-52] are shown. For the MRT curves the normalisation factors determined 
from the distributions are used. The FMS prediction is normalised to the fit result at Wr^p = 
90 GeV. 



26 



S100 



Q. 

—a 

T 



10 r 



1 r 



_ • H1 

; A ZEUS 

: }:^-=*^ 


(Q >, 
[GeV^] " 




a-^-rrJ"" 22.4: 


; <-=i^ 




: a) ■-- 


MRT (CTEQ6M) ; 


50 


100 300 

[GeV] 



1.5 - 



1 - 



0.5 - 



- 



• H1 

A ZEUS 



I 1 1 nil I I I I I I 



. 1 



b) 



0.1 



10 

[GeV^] 



Figure 6: a) Total cross sections for elastic J/tp production as a function of W-yp in the range 
\t\ < 1.2 GeV^ in electroproduction in three bins of Q^. (Q^) indicates the bin centre value in 
the range considered. The solid lines show fits to the HI data of the form a oc W^^^. The 
dashed curves show the MRT QCD prediction based on the gluon distribution CTEQ6M [48] 
with the normalisation factors from the fit to the distribution. Results from the ZEUS exper- 
iment [16] in a similar kinematic range are also shown. They have been scaled to the given {Q"^) 
values using the dependence measured by ZEUS, b) The fit parameter 5 as a function of Q^. 
The inner error bars show the statistical error, while the outer error bars show the statistical and 
systematic uncertainties added in quadrature. 
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Figure 7: Differential cross section dcr/dt for elastic J /ip production as a function of \t\ a) in 
four bins of Q"^ in the range 40 < W-^p < 160 GeV. {Q"^) indicates the bin centre value in the 
range considered. The inner error bars show the statistical error, while the outer error bars 
show the statistical and systematic uncertainties added in quadrature. The solid lines show fits 
to the data of the form da /At oc e^*. The dashed curve shows the result of a fit proposed by 
Frankfurt and Strikman [8]. Figures b) and c) show the photoproduction measurements in the 
ranges 135 < W^p < 235 GeV and 205 < W^p < 305 GeV. 
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Figure 8: Differential cross section da/dt for elastic J/^ production as a function of W^^ in 
bins of \t\ a) in photoproduction and b) in electroproduction. (Q^) and indicate the bin 
centre values in the ranges considered. The inner error bars show the statistical error, while the 
outer error bars show the statistical and systematic uncertainties added in quadrature. The solid 
lines show the results of one-dimensional fits of the form o oc W^t^^'^^-^^ in each t bin. 
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Figure 9: a) The effective trajectory a{t) as a function of \t\ in the range 40 < W^p < 305 GeV 
for photoproduction ((Q^) = 0.05 GeV^) and 40 < W^p < 160 GeV for electroproduction 
HQ"^) = 8-9 GeV^). The data points are the results of the one-dimensional fits shown in 
figure 8. The inner error bars show the statistical error, while the outer error bars show the 
statistical and systematic uncertainties added in quadrature. The solid and dashed lines show 
the results of two-dimensional fits (equation 2) together with la-error bands, which take the 
correlation between the fit parameters into account. A comparison with the results of the ZEUS 
collaboration [6, 16] is shown in b) and c) for photoproduction and electroproduction respec- 
tively. The data in [16] are derived at slightly different values of (Q^). The lines are results 
from the two-dimensional fits. 
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Figure 10: The values of the t slope parameter b{W^p) as a function of Wyp in the range \t\ < 
1.2 GeV^ for a) photoproduction and b) electroproduction. (Q^) indicates the bin centre value 
in the range considered. The data points are the results of one-dimensional fits of the form 
da/dt oc e^* in W^p bins. The inner error bars show the statistical errors, while the outer error 
bars show the statistical and systematic uncertainties added in quadrature. The solid lines show 
the results of the two-dimensional fits (equation 2) as in figure 9. In a) the data are compared 
with results from the ZEUS collaboration [6]. 
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Figure 11: The differential cross sections for diffractive (elastic and proton dissociative) J/tp 
production as functions of the four angles 9*, (p*, \1' and $ defined in the text. The data are 
shown in bins of Q'^ for the range 40 < W^p < 160 GeV and |t| < 5 GeV^. (Q^) indicates 
the bin centre value in the range considered. The inner error bars show the statistical error, 
while the outer error bars show the statistical and systematic uncertainties added in quadrature. 
The solid lines represent the results of fits to the data. The dashed lines in b) and d) are fits 
assuming s-channel helicity conservation. 
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Figure 12: Spin-density matrix elements as functions of (a-e) and \t\ (f-j) for the range 
40 < VKyp < 160 GeV. The data points are the results of fits of equations 4- 7 to the data shown 
in figure 1 1 . The inner error bars show the fit result including only the statistical error, while 
the outer error bars also include the systematic uncertainties. The expectations from SCHC are 
shown as solid lines. The results from the ZEUS collaboration are also shown, ( a), c) and 
f) [6,16] and b) [5,15]). 
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Figure 13: a) Ratio R = a^/a'^ as a function of Q"^ for the range 40 < Wyp < 160 GeV and 
It I < 5 GeV^. The data are compared with the result of a MRT calculation [26] based on the 
CTEQ6M [48] gluon distribution. Also shown are results from the ZEUS collaboration [6, 16]. 
b) The cross sections for longitudinally and transversely polarised photons and as a 
function of Q^. The MRT QCD calculations based on different gluon distributions ( [48-50]) 
are also shown with the same normalisation factors as derived from figure 4. The inner error 
bars show the statistical errors, while the outer error bars show the total errors. 
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Table 2: Cross section for the elastic process — > J/ipp measured in bins of for Wyp = 
90 GeV and for |t| < 1.2 GeV^. (Q^) indicates the bin centre value in the range considered. 
The first error is statistical and the second the total systematic uncertainty. 
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Table 3: Photoproduction cross section for the elastic process jp J/ipp in bins of W^p for 
\t\ < 1.2 GcV^ using the data sets II-IV (table 1). (Wyp) indicates the bin centre value in the 
W^p range considered. The first error on the cross section is statistical and the second the total 
systematic uncertainty. Note that there is an overlapping bin between data sets HI and IV at 
(W^p) = 219.6 GeV, which is averaged for figure 5 to cr = 151 ± 8 ± 20 nb. 
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[ GeV^] 


[ GeV^] 


[GeV] 


\ IP 1 

[GeV] 


a 

[nbl 


9 — 1^ 




zLD — 7n 


oo.o 


1 -i_ 9 Q -1- 9 4 






7n — inn 


oo.y 


Of) fiiq^ioQ 
ou.u zn o.^ zn z.c; 






100 - 130 


114.1 


41.5 ± 5.1 ± 4.0 






130 - 160 


144.2 


45.0 ±8.8 ±4.5 


5-10 


7.0 


40 - 70 


53.3 


12.9 ± 2.5 ± 1.2 






70 - 100 


83.9 


14.5 ±2.5 ±1.4 






100 - 130 


114.1 


24.7 ±4.1 ±2.4 






130 - 160 


144.2 


24.1 ±6.2 ±2.5 


10-80 


22.4 


40 - 70 


53.4 


3.19 ±0.69 ±0.31 






70 - 100 


83.9 


4.04 ±0.70 ±0.39 






100 - 130 


114.1 


5.29 ± 1.0 ±0.5 






130 - 160 


144.2 


6.10 ± 1.6 ±0.6 



Table 4: Total cross section for the elastic process jp J/ipp measured in bins of and W^p 
for \t\ < 1.2 GeV^. (Q^) and (W^p) are the bin centre values in the indicated ranges. The first 
error on the cross section is statistical and the second the total systematic uncertainty. 



g2[GeV2] (g2)[GeV^ 



6 [GeV 



-2 





<1 


0.05 


0.75 ±0.03 ±0.03 


4.57±0.06lg:i^ 


2 


- 5 


3.2 


0.67 ±0.20 ±0.14 


4.11 ±0.26 ±0.37 


5 


- 10 


7.0 


0.83 ±0.31 ±0.15 


3.50 ±0.50 ±0.49 


10 


- 80 


22.4 


0.69 ±0.32 ±0.14 


3.49 ±0.45 ±0.33 



Table 5: The parameters S (a oc W^p) and b oc e^*) measured in bins of in the range 
40 < W^fp < 160 GcV and \t\ < 1.2 GeV^. The values (Q^) indicate the bin centre value in the 
Q'^ range considered. The first error is statistical and the second systematic. 
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[GeV'] [GeV'] 


[GeV'] 


[GeV'] 


dcr/di 
[nb/GeV'] 


< 1 


0.05 





— 0.07 


0.03 


285 ± 9 ± 25 








— 014 


0.10 


180±7± 16 








— 91 


0.17 


130±6± 11 






n 91 


— 30 


0.25 


92.1 ±4.0 ±8.1 






0.30 


- 0.40 


0.35 


61.2 ±3.1 ±5.4 






n 4(1 


— fiO 


0.49 


32.5 ± 1.5 ± 2.9 






n fin 


— QO 


0.73 


10.60 ±0.60 ±0.90 






n Qo 


— 1 20 


1.03 


2.70 ±0.20 ±0.30 


2-5 


3.2 





- 0.08 


0.04 


107 ± 14 ± 10 






U.VJO 


— 018 

VJ. xo 


0.13 


95.1 ± 11.0 ±9.1 






18 


— 38 


0.27 


40.2 ±5.4 ±3.9 






38 


— 1 20 


0.68 


8.04 ± 1.05 ±0.77 


in 
— iU 


1 n 


U 


n nQ 
— U.Uo 


0.04 


78.6 ±13.2 ±7.5 






0.08 


- 0.18 


0.13 


27.7 ±5.7 ±2.7 






0.18 


- 0.38 


0.27 


18.9 ± 3.7 ± 1.8 






0.38 


- 1.20 


0.68 


5.21 ±0.96 ±0.50 


10-80 


22.4 





- 0.08 


0.04 


15.0 ±3.1 ± 1.4 






0.08 


- 0.18 


0.13 


8.90 ±2.14 ±0.85 






0.18 


- 0.38 


0.27 


4.55 ± 0.93 ± 0.44 






0.38 


- 1.20 


0.68 


1.36 ±0.25 ±0.13 



Table 6: Differential cross section for the elastic process 7p J /il)p measured in bins of 
and \t\ in the range 40 < W^p < 160 GeV using the data sets I and II (table 1). (Q^) and {\t\) 
are the bin centre values in the indicated ranges. The first error on the cross section is statistical 
and the second is the total systematic uncertainty. 





da/dt [nb/GeV^] 




|t| 0-0.07 


0.07 - 0.14 


0.14 - 0.30 


0.30 - 0.60 


0.60 - 1.20 


[GeV] 


[GeV2] 










45 


182 ± 20 ± 16 


115 ± 15 ± 10 


64.9 ±6.7 ±5.7 


35.5 ±3.6 ±3.1 


5.7 ±0.8 ±0.5 


55 


208 ± 20 ± 18 


118 ± 14 ± 10 


69.6 ±6.9 ±6.1 


35.6 ±3.4 ±3.1 


5.5 ±0.7 ±0.5 


65 


225 ± 23 ± 20 


169 ± 18 ± 15 


107.1 ±9.4 ±9.4 


34.3 ±3.7 ±3.0 


6.2 ±0.9 ±0.5 


75 


321 ± 31 ± 28 


151 ±19 ±13 


93.4 ±9.5 ±8.2 


38.9 ±4.3 ±3.4 


7.1 ± 1.1 ±0.6 


85 


292 ± 29 ± 26 


178 ± 20 ± 16 


132 ± 11 ± 12 


41.4 ±4.4 ±3.6 


8.2 ±1.1 ±0.7 


95 


326 ± 32 ± 29 


224 ± 24 ± 20 


135 ±12 ±12 


46.3 ±4.8 ±4.1 


7.5 ±1.1 ±0.7 


105 


392 ± 37 ± 34 


224 ± 26 ± 20 


125 ± 12 ± 11 


48.5 ±5.3 ±4.3 


7.7 ±1.1 ±0.7 


119 


376 ± 31 ± 33 


265 ± 24 ± 23 


142 ± 11 ± 12 


60.9 ±4.9 ±5.4 


8.2 ± 1.0 ±0.7 


144 


458 ± 42 ± 40 


267 ± 29 ± 23 


167 ± 14 ± 15 


51.5 ±5.3 ±4.5 


8.4 ± 1.1 ±0.7 


181 


537 ± 28 ± 68 


427 ± 18 ± 46 


202 ± 9 ± 25 


67.0 ±4.1 ±8.8 


10.7 ± 1.1 ± 1.7 


251 


744 ± 48 ± 88 


573 ± 28 ± 75 


246 ± 13 ± 32 


71.6 ±5.5 ±9.8 


13.1 ± 1.8 ± 1.9 



Table 7: Differential photoproduction cross sections da/dt for the elastic process jp — > J/ipp 
measured in bins of W^p and \t\ using data sets II-IV (table 1). The first error is statistical and 
the second the total systematic uncertainty. 
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da/dt [nb/GeV^ 




(W^,)[GoV] 




t 


[ GeV-] 0-0.1 0.1-0.3 0.3-1.2 


57 
98 
140 


33.3 ±4.9 ±3.2 17.3 ±2.2 ±1.7 3.4 ± 0.5 ± 0.3 
51.3 ±6.7 ±4.9 30.1 ±3.5 ±2.9 5.5 ±0.7 ±0.5 
60 ±12 ±6 31.5 ±5.8 ±3.0 6.0 ± 1.1 ± 0.6 



Table 8: Differential electroproduction ({Q^) = 8.9 GcV^) cross section da/dt for the elastic 
process J/ipp measured in bins of W^p and \t\ using data set I (table 1). The first error is 

statistical and the second the total systematic uncertainty. 



[GeV'] 


\t\ 

[GeV'] 


(l*l> 

[GeV'] 




0.05 


- 0.07 


0.03 


1.202 ±0.012 ±0.017 




0.07-0.14 


0.10 


1.240 ±0.012 ±0.019 




0.14-0.30 


0.22 


1.195 ±0.011 ±0.016 




0.30 - 0.60 


0.43 


1.121 ±0.013 ±0.018 




0.(30 - 1.20 


0.81 


1.117 ±0.021 ±0.019 


8.9 


- 0.1 


0.05 


1.173 ±0.064 ±0.038 



0.1 - 0.3 0.19 1.185 ±0.054 ±0.037 
0.3 - 1.2 0.64 1.168 ±0.059 ±0.037 

Table 9: The effective Pomeron trajectories a{t) derived from one-dimensional fits of the W^p 
dependence in bins of t, in the ranges Q' < 1 GcV', 40 < W^p < 305 GcV (photoproduction) 
and 2 < g' < 80 GeV', 40 < W^p < 160 GeV (electroproduction). The values (Q') and {\t\) 
are the bin centre values in the indicated ranges. The first error is statistical and the second 
systematic. 
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Data set 


[GeV^] 


[GeV] 


[GeV] 


[GeV"'] 


I 


8.9 


40 - 80 


57.3 


3.77 ±0.34 ±0.33 






80 — 120 


98.2 


3.79 ± 0.29 ± 0.32 








lOc/.U 


o.o^ HZ u.^o m u.oo 


II 




40 — 50 


44 8 


4 13 ± 20+oJ,^ 






50 — 60 

t_/ W \J\J 


54 8 

'.J t: . 


4 30 ± 19+S-if 






60 - 70 


64.8 


4.57±0.20+o}7 






70 — 80 


74 8 


4 46 ± 24+S i5 

^ ^ ^ — 0.46 






cf) — 90 


84 9 


^■^'J ^ '-'■■^'-'-0.25 






Qn — 1 nn 


ij'-±..ij 


4 79 _|_ n 91 +015 






1 nn 1 1 n 

lUU — IIU 


1 OA Q 


4 7Q -1- n 99+0.15 






110 - 130 


119.5 


4.71±0.16i|J;i^ 






130 - 160 


111.1 


1.95 ±0.19l[5i^ 


ni 


0.05 


135 - 235 


180.6 


5.08±0.14l°|^ 


IV 


0.05 


205 - 305 


250.7 


5.41 ± 0.20l^i^ 



Table 10: The slope parameter h derived from one-dimensional fits to the t dependence mea- 
sured in bins of W^p. The values (Q^) and (VF-yp) are the bin centre values in the indicated 
ranges. The first error is statistical and the second systematic. 



m 

[GeV^] 


^04 
'00 


T.04 
'1-1 


R 









0.05 -0.030 ±0.016 ±0.027 0.020 ± 0.016 ± 0.042 -0.030+°;°}^!°;°^^ 

3.2 0.049 ± 0.079 ± 0.050 -0.129 ± 0.070 ± 0.039 0.052;°;°^™! 

7.0 0.19 ±0.14 ±0.06 -0.017 ±0.10 ±0.04 0.23l[!:?^l[!:J^^ 

22.4 0.38 ±0.16 ±0.06 -0.04 ± 0.12 ± 0.04 0-62+° j^l°;}^ 

Table 11: The spin-density matrix elements, Tqq and r^i^, and the ratio of cross sections of 
longitudinally and transversely polarised photons i? as a function of in the range \t\ < 
5 GeV^ and 40 < W^p < 160 GeV. The values (Q^) indicate the bin centre values in the 
range considered. The first error is statistical and the second systematic. 



[GeV^] 


rU 


^00 + 2rJ]^ 


7-5 1 9„5 
'00 ^ ^' 11 


3.2 


0.149 ±0.077 ±0.064 


0.026 ±0.035 ±0.026 


-0.035 ±0.072 ±0.055 


7.0 


0.43 ±0.11 ±0.06 


0.062 ±0.054 ±0.028 


-0.16 ±0.12 ±0.06 


22.4 


0.53 ±0.13 ±0.05 


0.026 ±0.069 ±0.031 


0.04 ±0.16 ±0.06 



Table 12: The spin-density matrix element rj^^^ and the combined elements Tqq + 2r|^ and 
r^Q + 2rfi as a function of in the range |t| < 5 GeV^ and 40 < W^p < 160 GeV. The values 
(Q^) indicate the bin centre value in the Q'^ range considered. The first error is statistical and 
the second systematic. 
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(Q^) [GeV^ 


{\t\)[GeV' 


T.04 
'00 


^.04 
'1-1 


0.05 


0.03 


0.003 ± 0.039 ± 0.028 


-0.011 ±0.036 ±0.030 




0.10 


0.011 ±0.043 ±0.029 


-0.041 ±0.042 ±0.030 




0.22 


0.026 ±0.036 ±0.028 


0.104 ±0.035 ±0.029 




0.43 


0.013 ±0.037 ±0.029 


0.025 ±0.037 ±0.030 




0.84 


0.047 ±0.041 ±0.029 


0.064 ± 0.047 ± 0.034 




1.8 


0.066 ±0.061 ±0.028 


-0.010 ±0.060 ±0.030 




3.5 


0.018 ±0.081 ±0.028 


-0.074 ±0.082 ±0.032 


8.9 


0.05 


0.32 ±0.15 ±0.06 


-0.13 ±0.11 ±0.04 




0.19 


0.22 ±0.13 ±0.06 


-0.07 ±0.10 ±0.04 




0.64 


0.05 ±0.10 ±0.05 


-0.071 ±0.083 ±0.036 




3.0 


0.23 ±0.19 ±0.06 


0.06 ±0.13 ±0.03 



Table 13: The spin-density matrix elements, Tqo and r^ti, as a function of \t\ in the range 
40 < W^p < 160 GeV for photoproduction and electroproduction. {Q"^) and {\t\) are the bin 
centre values. The first error is statistical and the second systematic. 



( 


t 


) 


[GeV^ 


-1-1 


-00 + 2rli 


-00 + 2r5, 



0.05 0.19 ±0.11 ±0.06 -0.16 ±0.12 ±0.06 0.031 ± 0.055 ± 0.029 

0.19 0.62 ±0.11 ±0.06 0.04 ±0.10 ±0.06 0.004 ± 0.048 ± 0.026 

0.64 0.361 ±0.097 ±0.061 -0.073 ± 0.091 ± 0.057 -0.039 ± 0.043 ± 0.027 

3.0 0.07 ±0.15 ±0.06 -0.15 ± 0.15 ± 0.06 0.083 ± 0.079 ± 0.030 



Table 14: The spin-density matrix element rl_i and the combined elements Tqq ± 2rli and 
r^o ± 2rfi as a function of \t\ in the range 40 < W^p < 160 GeV and 2 < < gg GeVl 
indicates the bin centre value. The first error is statistical and the second systematic. 
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